Collective behavior in a granular jet: Emergence of a liquid with zero surface-tension 
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We perform the analog to the "water bell" experiment using non-cohesive granular material. 
When a jet of granular material, many particles wide, rebounds from a fixed cylindrical target, it 
deforms into a sharply-defined sheet or cone with a shape that mimics a liquid with zero surface 
tension. The particulate nature of granular material becomes apparent when the number of particles 
in the cross-section of the jet is decreased and the emerging sheets and cones broaden and gradually 
disintegrate into a broad spray. This experiment has its counterpart in the behavior of the quark- 
gluon plasma generated by collisions of gold ions at the Relativistic Heavy Ion Collider. There a 
high density of inter-particle collisions gives rise to collective behavior that has also been described 
as a liquid. 

PACS numbers: 45.70.-n, 45.70.Mg, 47.50.Tn 
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When one or two particles strike a smooth wall at nor- 
mal incidence, they rebound in the direction whence they 
came. Yet, as we show here, a dense stream of non- 
cohesive particles hitting a target retains its integrity and 
deforms into a thin sheet with a shape resembling the 
structures created by an impinging water jet 1, 2]. Thus 
the collective behavior of many particles differs qualita- 
tively from that of the individual components. When 
can a jet of discrete particles be modeled as a liquid 
[1) 0) S @| and how do the liquid patterns emerge out of 
individual particle scattering events from a target? Such 
questions, posed and investigated here with granular ma- 
terials, have their counterpart in much more microscopic 
situations such as the quark-gluon plasma caused by rel- 
ativistic high-energy collisions of gold ions, which also 
produces scattering patterns indicative of a liquid state 
(author?) 0, i]. Our findings provide a macroscopic, 
purely classical example of how strong interactions, me- 
diated by rapid collisions in a densely packed region, can 
give rise to liquid behavior. 

When a water stream hits a flat target, it spreads sym- 
metrically in the direction transverse to the impact and 
deforms into a thin sheet. For targets smaller than the 
stream diameter, the sheet forms a hollow bell-shaped 
structure that envelopes the target. Such "water bells" 
were reported in 1883 by Savart [l| and have been sys- 
tematically studied by Clanet 01 ■ By contrast, with the 
exception of several simulations and experiments on gran- 
ular streams passing by solid obstacles at relatively low 
particle density and speed [13, Hi, little is known 
about the structures emerging from equivalent experi- 
ments using granular material rather than a liquid. Here, 
we examine a non-cohesive granular jet colliding with a 
target. For jets many particles wide, the impact pro- 
duces granular sheets and cones similar to those seen for 
water with structures that depend on the ratio of the jet 
to the target diameter. In this regime, we conclude that 
the granular medium behaves as if it were a liquid with 
infinite Weber number - that is, with zero surface ten- 



sion appropriate for a system of non-cohesive particles. 
The particulate nature of the material becomes appar- 
ent when the sheets and cones broaden and gradually 
disintegrate as the number of particles in the beam is 
decreased. 

We prepared our dense jets by packing granular ma- 
terial into a 40cm section of a glass launching tube of 
inner diameter Dj^t = 0.73cm (FiglT^). The grains, 
which were compacted to a reproducible density by tap- 
ping, were mono-disperse spherical beads of glass (p = 
2.5g/cm'^) or copper (p = 8.2g/cm'^) with diameters be- 
tween d — 50/im and d = 2.1mm. Prior to filling the 
tube, we baked the beads in vacuum to minimize any 
residual adhesion between beads. Pressurized gas accel- 
erated this granular plug into a jet which hits a target 
2.5cm in front of, and collinear with, the tube. The jet ve- 
locity Uq could be varied between 1 and 16m/s. The im- 
pacts were filmed at 2000 frames/second with a Phantom 
V7.1 camera. Fig[T|3 and c show the side and front views 
respectively of a colliding jet. Since there are very few 
collisions inside it, the jet maintains its cylindrical shape 
before hitting the target, consistent with jets produced 
when a sphere impacts loosely packed powers [12, [l3| . 
After hitting the target, the jet deforms into an extraor- 
dinarily thin symmetric granular sheet clearly resembling 
a spreading liquid. To investigate the similarities be- 
tween granular and ordinary liquid jets, we first keep the 
jet diameter, Dj^t, fixed at values much larger than the 
particle diameter, d, and vary the target diameter, Dxan 
as shown in Fig[TJi. When Dxar is reduced, these planar 
sheets change into cones with opening angle, Vo < 90°. 
FigE] shows that ipo increases linearly with Dxar/ Djet 
until ipo saturates at 90° above Dxar/ Djet ~ 2. Glass 
and copper beads produces essentially identical behavior. 

We can compare our results with those from Clanet 
for water jets, who found that the opening angle, 
■0, of the "water bell" depends on the Weber number 
We — pU§Djet/cr, where Uq is the jet velocity and p 
the density and a the surface tension of water. In the 
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FIG. 1: (color online). Liquid-like granular sheets and cones, 
(a) Sketch of the experimental setup, (b) Side view of a gran- 
ular jet comprised of 100/im glass beads hitting a target at 
velocity Uo ^ lOm/s. The ratio of target to jet diameter is 
Drar / Djet ~ 3.80. The images show the jet 0.5ms before as 
well as 2.5ms and 9.5ms after impact (left to right). Scale bar 
is 3.0cm. (c) Axial view of the same granular sheet 35.5ms 
after impact. Scale bar is 2.5cm. (d) Side views of the hollow 
cones produced by lOO/xm glass bead jets for smaller target- 
to-jet diameter ratios. Drar/ Djet = 2.0, 1.62, and 0.88 (left 
to right). Scale bar is 6.0cm. 



large- We limit (large ?7o or small ct), -0 approaches a 
constant, Vo, that depends only on DTar/ Dj^t- Fig[2] 
shows excellent overlap of our values for ipQ with those 
of Clanet. As shown in the inset to Figl^l ij] for granular 
material remains constant as L/q is varied over our entire 
experimental range (an order of magnitude) whereas V' 
for water decreases at small velocities. This suggests a 
granular jet impacting a the target behaves like a liquid 
with negligible surface-tension: regardless of the value of 
L/o, the Weber number is pinned at infinity since the ab- 
sence of significant cohesive forces effectively drives the 
surface tension to zero. 

The linear dependence of ipQ on Dxar / Djet in the 
large- We limit can be understood from momentum con- 
servation. The magnitude of the momentum of the in- 
coming granular jet reaching the target in time t is 
P^^ = (7r/4)Djgj[/Qr. The collisions of the jet with the 
target are inelastic so the outgoing momentum magni- 
tude Pout = CPin where C is the average coefficient of 
restitution. Along the axis, the momentum balance is 

given by: P,n ~ Pout COSV'O = -Pm(l - CCOS V'o) = FrarT, 
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FIG. 2: Opening angle of the granular sheets and cones. An- 
gle t/iQ as a function of Dxar / Dj^t for lOO/^m glass beads (■), 
100/im copper beads (v) and water in the large-We limit (o). 
The only observable difference is that tpo at large Drar / Djet 
is ~ 85° for water but 90° for grains. The solid lines are fits 
of Eq[T]with {A ~ B) ^ 0.30 for DTar/Dj^t < 2. Inset: An- 
gle -0 as a function of jet velocity Uo for 100/im glass beads 
with DTar/Dj,t = 2.47 (■) and Drar/Dj^t = 1-32 (.), and 
for water with Drar / Djet = 3.29 (a) and Dxar / Djet ~ 1-33 
(0). The horizontal lines are %p = 90° (top) and tp = 65.5° 
(bottom). Data for water is taken from 



where Fxar is the average force along the axis exerted 
by the target on the jet. Transverse to the axis, the mo- 
mentum (averaged over the entire sheet) remains zero. 
When DTar/Dj,t < 1, ^Va. = {ATr/A)D^^^pUi. Here A 
is a constant describing the average glancing collision an- 
gle for a particle. In this region, the average restitution 
coefficient, C, also depends on the fraction of particles 
hitting the target: C = [1 - B{DTar/ Djetf] where B 
depends on the coefficient of restitution for single colli- 
sions. When Dxar/ Djet S> 1, the entire momentum of 
the jet is reflected by the target. Finally we obtain: 
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Clanet reached the same result by considering the mo- 
mentum transfer during the impact of a water jet on a 
target using hydrodynamic equations 0. By fitting the 
experimental data in FigH we find (A-B) = 0.30 ± 0.02 
which is close to the value « 0.352 for water 

The impact of dense particle streams clearly generates 
similar patterns as do liquids. How does the particu- 
late nature of granular material becomes manifest as the 
number of particles within the jet decreases? To vary this 
number we change the ratio, Djet/d, of jet to particle di- 
ameter. In the limit DTar/ Djet ^ 1, FigsO^-c show 
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major qualitative changes in the particle trajectories as 
Djet/d is reduced. The images were created by superim- 
position of many different consecutive still images; each 
pixel shows the maximum intensity at that location over 
all images in the time period. For Djet/d = 73, almost 
all particles emerge in a sheet normal to the jet axis. For 
Djet/d = 14.6, the sheet becomes more diffuse as some 
particles leave the plane as shown by the bright lines. For 
Djet/d — 3.5, a firework-like pattern results after impact 
and sheet structure is no longer apparent as particles re- 
bound from the target in a broad angular distribution. 
To quantify this trend, we plot in Fig[3li the angular 
scattering distribution, obtained by averaging images as 
in FiglSK-c along a circle of radius r = 8.7 D jet centered 
on the jet axis (shown partially in FiglS^a) as the dashed 
line). The axial position of the center is chosen to be 
in the middle of the ejected particles, which gradually 
moves away from the target in upstream direction as the 
particle diameter d increases. To avoid the target holder, 
only the scattering profile along the upper half of the 
circle is plotted. Each profile is normalized to its peak 
value. The scattering angle 9 is zero along the axis of the 
jet and increases clockwise. As Djet/d decreases, these 
profiles become broader. 

To investigate this crossover to diffuse scattering we 
examined the roles of air and inelasticity during inter- 
particle collisions. To exclude air as the primary cause of 
the sheet formation, three experiments were performed, 
(i) Doing the entire experiment in a helium atmosphere 
and using helium as the accelerating gas produced no 
change in the angular distribution, as shown in Figl3jl 
for Djet/d — 73. (ii) Performing our experiment at re- 
duced ambient pressure of 31Pa also resulted in a thin, 
planar sheet structure. (31Pa corresponds to a mean 
free path of air molecules three times the 100/zm grain 
diameter used in FiglHK. Although air accelerates the 
granular column from behind, we estimate that the ve- 
locity of air penetrating into a 40cm granular pack com- 
prised of lOO/im beads with packing fraction around 0.6 
is two orders of magnitude smaller than the velocity of 
the jet. Thus, by the time the front of the jet hits the 
target, a negligible amount of air has entered the cham- 
ber.) (iii) We decreased the tube diameter Djet, but 
fixed the particle-air interaction by keeping d = 100/im. 
In this case the planar structure shown in Figl3^ dis- 
appeared gradually, and for Dj^t ^ d we regained the 
firework pattern shown in Figl3l:. In addition, we altered 
the inter-particle collision dynamics by changing the par- 
ticle material or the surface roughness. Fig[3Ji shows that 
for copper particles that are less elastic than glass or for 
rough sand particles the angular scattering distribution 
remains unchanged. However, distributions for the same 
Djet/d do not overlap quantitatively for different d, and 
the liquid-like sheets tend to be more sharply delineated 
with smaller particles. 

While we cannot rule out that inelastic collisions or 




FIG. 3: (color online). Scattering patterns as a function 
of particle density across the jet diameter, (a)-(c) Maxi- 
mum intensity projections of the scattering process for fixed 
Drar / Djst = 3.8 but varying particle density: (a) Djet/d = 
73; (b) Djet/d = 14.6; (c) Djet/d = 3.5. (d) Normalized scat- 
tering profiles for glass beads in air with Dj^t/d = 73 (■), 
Djet/d = 29.2 (o), Dj,t/d = 14.6 (A), Djet/d = 3.5 (V); 
glass beads in helium Djet/d — 73 (red ♦); copper beads in 
air Djet/d = 73 (green <) and rough sand in air Djet/d = 73 
(blue 0). The arrow and the dashed circle in the inset indicate 
the direction and the cross-section of the jet. (e) Half-width 
^haif of the angular scattering profiles in (d) as a function of 
d/Djet- Solid line is a linear fitting. The model is sketched 
in the inset and explained in the text. 



effects of air affect the detailed shape of the scatter- 
ing profiles, our results imply that the creation of the 
sheets arises fundamentally from the rapid collisions oc- 
curring in an interaction region right in front of the tar- 
get. We found similar behavior when two granular jets 
collide head-on, implying that the target serves primarily 
to reverse the direction of particles incident upon it. For 
this situation, a simple geometric model can capture the 
essence of the crossover from fluid- to granular behav- 
ior. As shown in FiglSJj, inset, we divide the system into 
three zones: two external regions in which the jets are 
traveling towards one another but have not yet collided, 
and an interaction region with the same diameter as the 
jets, -Djet, and axial length, /. Inside the interaction re- 
gion particles undergo rapid collisions and are confined 
by pressure from the incoming jets on both sides. The 
only way in which particles can escape (which they must 
since new particles are entering the region continuously) 
is to emerge perpendicularly to the jet axis. A measure 
for deviations from this transverse axis and thus for the 
half- width of the angular distributions in FigjSli is the an- 
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FIG. 4: Azimuthal scattering profiles for granular jets with 
rectangular cross-section (aspect ratio 2). Shown are views 
along the axial direction, looking at the target from behind, 
(a) Anisotropic profile resulting for Dj^t/d = 73, where Dj^t 
is the long axis of the jet cross-section. The target holder 
blocks the view of the scattered portion shooting downward. 
Scale bar is 3.0cm. (b) Isotropic profile resulting for Djet jd — 
7. (c) Azimuthal profiles along a circle with r = 8.7Djet for 
Djet/d = 73 (■), Dj,t/d = 29.2 (o), Dj,t/d = 14.6 (a), 
Dj^tjd = 7 (V) and Dj^t/d = 3.5 (♦).The direction and the 
cross section of jet are shown in the inset; is taken to be 
zero along the long axis of the jet cross-section. 



patterns for granular jets with rectangular cross-section. 
To quantify the anisotropy and compare with the RHIC 
results 0, BIlll' analyze the second coefficient, V2, of 
the Fourier expansion of the azimuthal scattering profiles 
in FigO:. For Djet/d = 73 we obtain 1/2 — 0.16 which is 
larger than the typical value at the RHIC experiments, 
indicating a stronger interaction between component par- 
ticles. As FigsHb and c show, the scattering patterns 
become more isotropic with decreasing Dj^t / d. 

These experiments demonstrate how non-cohesive 
granular material can produce collective motion reminis- 
cent of liquids. Our results suggest that rapid particle 
collisions in a very narrow interaction region, with thick- 
ness of a few particle diameters, eject very thin granular 
sheets that mimic the liquids without surface tension. 
The crossover between diffuse and sharp scattering pro- 
files appears to be controlled primarily by the number of 
collisions in the jet. These findings may have analogs in 
other, disparate parts of physics, where a high density of 
collisions dominates the behavior. 

We thank R. Bellwied, E. Corwin, S. Gavin, T. 
Poschel, J. Royer, T. Witten and L. Xu. This work was 
supported by the NSF through its MRSEC program un- 
der DMR-0213745 and through its Inter- American Ma- 
terials Collaboration under DMR-0303072. 



gle 0haif by which particles can escape ballistically from 
the center of the beam. By geometry, this is given by 
tan(6'haif) = l/{Djet)- For dense particle streams where 
the mean free path is much smaller than the diameter of 
the particles, we expect / ~ 2d since that is the smallest 
region occupied by the two colliding particles. Therefore, 
the half-width should scale as ^haif = arctan(A2(i/£)jet) 
where A is a constant 0(1). This dependence on parti- 
cle density d/Djet across the beam is indeed born out 
by the data in Figllh, with A = 1.05 ± 0.05. As a result, 
when Djet /d is sufficiently large the majority of particles 
undergo rapid, multiple collisions in front of the target 
and are ejected into a narrow angular, similar to a liquid 
film. For example, 6'haif < 2° for jets with Djet/d > 60. 

When the granular scattering mimics that of a liquid, 
we expect azimuthally anisotropic patterns to be created 
by jets with rectangular cross-sections. The incoming jets 
establish an elongated interaction region within the plane 
transverse to the beam (azimuthal 0-direction) so that 
more particles will escape along the direction of the short 
axis of the cross-section where a greater pressure gradi- 
ent exists. This is also the situation pertinent to non- 
central collisions of heavy ions at the Relativistic Heavy 
Ion Collider (RHIC) where partial overlap of the collid- 
ing ions establishes an elongated interaction zone for the 
quark-gluon plasma. Azimuthally anisotropic scattering 
patterns have been taken as the main evidence for liquid- 
like behavior of this plasma 0, @]- As shown in Fig [3^ 
for Djet/d ^ 1, we observe sharply focused azimuthal 
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